ON A PROBLEM OF ERDOS INVOLVING 
THE LARGEST PRIME FACTOR OF n 

Aleksandar Ivic 

Abstract. Let P(n) denote the largest prime factor of an integer n (> 2), and let 
N{x) denote the number of natural numbers n such that 2 < n < a;, and n does not 
divide P{n)\. We prove that 

N{x)=x(^ + o(^^\og2xl\ogxy^ p{\ogx/\ogt)^^dt, 

where p{u) is the Dickman-de Bruijn function. In terms of elementary functions we 
have 

N{x) = a; exp I - >/2 log a; log2 a; (1 + ©(logg x/ logg a;)) | , 

thereby sharpening and correcting recent results of K. Ford and J.-M. De Koninck 
and N. Doyon. 



1. Introduction and statement of results 

Let P{n) denote the largest prime factor of an integer n (> 2). In 1991 P. Erdos 
[5] proposed the foUowing problem: prove that N{x) = o{x) (x oo), where N{x) 
denotes the number of natural numbers n such that 2 < n < x, and n does not 
divide P{n)\. This problem is connected to the so-called Smarandache function 
S{n), the smallest integer k such that n\k\. 

Erdos's assertion was shown to be true by I. Kastanas [13], and S. Akbik [1] 
proved later that N{x) <^ xexp(— i A/log x ) holds. K. Ford [8] proposed an as- 
ymptotic formula for N[x). His Theorem 1 states that 

(1.1) N{x) ~ ^/^li±i^(logxlog2x)3/4a;i-i/-VK) (x oo). 

In this formula one should correct the constant, as will be show in Section 4. Here, 
as usual, logj^ x = logx and log^^ x = log(log;j_j^ x) for k > 2. The function p{u) is 
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the continuous solution to the difference delay equation up'{u) = —p{u — 1) with 
the initial condition p{u) = 1 for < tt < 1 and is commonly called the Dickman 
(or Dickman-de Bruijn) function. The function uq = uo{x) is implicitly defined by 
the equation 

log a; = Uq (^x^^'^o - , 

so that 

(1.2) uoix) = ( ^y^' 1 1 _ + + o ( 



^og2X J y 21og2a; 21og2a: \^Vlog2a^ 

The asymptotic formula (1.1) and the explicit expression (1.2) enabled Ford (op. 
cit.) in Corollary 2 to deduce that 

(1.3) N{x) = xexp |-(V2 + o(l))A/loga;log2a;| {x ^ oo). 

Unaware of Ford's work, J.-M. De Koninck and N. Doyon [3] recently tackled 
this problem again. In [3] they published the result that 

(1.4) N{x) — xexp | — (2 + o(l))-^loga;log2 x| (x — >• oo). 

Unfortunately, (1.4) is not true (it contradicts (1.3)). Namely the argument in 
[3] leading to their (3.8) does not hold, and the crucial parameter u = log x/ logy 
in the formula for \l/(a:,y) has to be evaluated more carefully. A straightforward 
modification of their proof leads then again to (1.3). 

The aim of this note is to sharpen (1.1) and (1.3), and to provide asymptotic 
formulas for sums of S'^{n) and 1/5''" (n) when r e R is fixed. The results are 
contained in the following 

THEOREM 1. 

(1.5) N{x) = xe■x.^\^-^/2L{x) |^l + ^o(a^) + 
where 



(1.6) L = L{x) = -\/log X log log X, 

and, for r > —1, 

_ log3a: + log(l + r)-2-log2 / 
(1-7) ^ , , 



2 log2 X V logs ^ 

(log3 X + log(l + r) - log 2) 



8 log2 X 



On a problem of Erdos 



3 



THEOREM 2. 



(1.8) N{x)=x\2 + 0{,ri^ W I Pl'^l^dt. 





2 "Viogt; 



THEOREM 3. Let r > be a fixed number. Then we have 
(1-9) Yl ^^ = xexp\-V2^L{x){l + gr-i{x) + 

2<n<x ^ ^""^ [ \ 

where L{x) and gr{x) are given by (1.6) and (1.7), respectively. Moreover, for 
any fixed integer J > 1 there exist computable constants ai^r = Ci''^ + 1)? 
a2,r • • • 5 CLj,r suck that 

(1.10) y s'in) = x'-+^ y +o( -^^^] . 

2<n<x 3=1 ° 

The asymptotic formula (1.8) is sharper than (1.5)-(1.7), since it is a true as- 
ymptotic formula, while (1.5)-(1.7) is actually an asymptotic formula for log ^^^^ . 
On the other hand, the right-hand side of (1.5) is given in terms of elementary 
functions, while (1.8) contains the (non-elementary) function p{u). Thus is seemed 
in place to give both (1.5) and (1.8), especially since the proof of (1.5) requires 
only the more elementary analysis which follows the work of De Koninck-Doyon 
[3]. 

The asymptotic formula (1.9) seems to be new, while (1-10) sharpens the results 
of S. Finch [7], who obtained (1.10) in the case J = 1. 

Acknowledgement. I thank Prof. G. Tenenbaum for pointing out formula 
(2.6) to me, and indicating how it improves the error term 0(log3 x/\og2 x) in the 
first version of Th. 2 (coming from (2.5)), to O ( \/\og2 x / log x) . I also thank the 
referee for valuable remarks. 

2. The necessary results of ^'(a;,y) 
The proofs of both Theorem 1 and Theorem 2 depend on results on the function 

(2.1) *(a;,y) = J] 1 (2 < y < x), 

n<x,P{n)<y 
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which represents the number of n not exceeding a;, all of whose prime factors do 
not exceed y. All the results that follow can be found e.g., in [9], [10] and [14]. 
We have the elementary bound 

(2.2) ^{x,y) < xexp ( - j"^^ ) (2 < y < x) 



2 logy 

and the more refined asymptotic formula 

(2.3) ^ ^' '\ V logy 

u — logx/logy, exp((log2 x)^/"^"''^) <y <x. 

Note that the Dickman-de Bruijn function p{u) admits an asymptotic expansion, 
as M — > oo, which in a simplified form reads 




p{u) — exp ^logu + log2 u — 1 + O 
Let, for given w > 1, ^ = ^(u) be defined by = - 

(2.4) = log « + log2 u + + O 

If u > 2, \v\ < ^u, then we have the asymptotic formula 

(2.5) p{u-v) = p{u)e-x^{vi{u) + 0{{l + v'^)/u)] . 

This result (see e.g., Hildebrand-Tenenbaum [10]) for bounded \v\ (see (3.3.9) of 
de la Breteche-Tenenbaum [2] ) yields the asymptotic formula 

(2.6) p(«-^;) =p(«)e<(") • |l + (^^^1 (|^| < ^o)- 

3. Proof of Theorem 1 

We turn now first to the proof of Theorem 1. The idea is to obtain upper 
and lower bounds of the form given by (1.5)-(1.7). For the former, note that if 
P^(n)|n, then n is counted by N{x), since P^{n) cannot divide P(n)!. Therefore, 
if Tq{x) denotes the number of natural numbers n < x such that P^(n)|n, then 
N{x) > Tq{x). The desired lower bound follows then from the formula for To{x), 
which is a special case of the formula for (r > is a given constant) 

(3.1) Tr{x) := Yl 



2<n<a;,p2(n)|n ^ ' 
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obtained by C. Pomerance and the author [12], and sharpened by the author in 
[11]. Namely for Tr{x) it was shown that 

(3.2) Tr{x) = xexp |-(2r + 2)^/^L{x) + gr{x) + O 

where L{x) and gr{x) are given by (1.6) and (1.7). We note that Erdos, Pomerance 
and the author in [6] investigated the related problem of the sum of reciprocals of 
P{n). They proved that (this result is also sharpened in [11]) 

The same argument leads without difficulty to 

,3.4) r,,.) . (^4,.) ^4^0 [^)) j\ (^) if 

It remains yet to deal with the upper bound for N{x). As in [3], we shall use 
the inequality {p denotes primes) 

(3.5) N{x)< E 

2<r<logK/log2 p<xi/r- ^ 

To see that (3.5) holds note that, if n is counted by N{x), then n must be divis- 
ible by which does not divide P{n)l. The condition r > 2 is necessary, since 
squarefree numbers q (> 1) divide P(g)!. Moreover, pr < P{n) cannot hold, since 
otherwise the numbers p, 2p, ... ,rp would all divide P{n)\, and so would p^, which 
is contrary to our assumption. Thus n — p^m (r > 2), P{m) < P{n) < pr, which 
easily gives (3.5). 

We restrict now the ranges of r and p on the right-hand of (3.5) to the ones 
which will yield the largest contribution. 
The contribution of r > 3L is trivially 

•C ^ y^xp""" < x2~^^ = a;exp(-31og2 • L), 

r>3L p 

which is negligible, since 3 log 2 > ^/2. Similarly, the contribution of p > exp(2L) 
is negligible, and also for p > exp((logx)^/'^) and r > (logx)^/^ log2 a; the contri- 
bution is, for some C > 0, 

<^ xp~'^ <C a;exp(— C\/loga; • log2a;), 

r>(log x)^/^ log2 X j?>exp((log x)^/^) 
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which is neghgible. 

With (3.5) and (2.3) we see that in the range 



Li :— exp((loga;)^/"^) <p< exp(iL), r < (loga;)^/^ log2 x 

we have 

log(.r///) log, 7- - rlogp 

^ — ~i — ~i 1 

log pr log p + log r 

^ log. + 0((logx)V3+.) ^ log. . ((log,,.-i/3)) . 
logp + 0(log2a;) logp V \\ ^ ^ jj 

Therefore by (2.3) we obtain, as a; — > oo, 

log X f log X 



^ — ,pr < — exp - l + o l))-^log 
\p^ ) p \ logp 



logp 



which gives 

' X 



Li<p<exp(iL) r>2 



p^ 



pr 

<.x > max exp -(l + o(l))- log' 



Li<p<exp(iL) 

< a;exp(-(2 + o(l))L), 



logp \logp 



which is negligible. The contribution of p < Li, r > (loga;)^/^ log2 x is easily seen 
to be also negligible. This means that the main contribution to the right-hand 
side of (3.5) comes from the range 

(3.6) exp(iL) < p < exp(2L), 2 < r < (logx)^/^ loga 

or more precisely, for some C > -\/2, 
N{x) < xexp(-CL) 

(3.7) + J2 E ^(§^^p^ 

2<r<(loga;)V6 log^ x exp(|L)<p<exp(2L) 

By trivial estimation it transpires that the contribution of r > 8 in (3.7) will be 
negligible, so that we may write 

(3.8) N{x) < X exp(-CL) + L max y Tc^rip^)-! 

c>0 ' 

2<r<7 
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where, for a given constant c > and r e N, 



exp(cL— l)<p<exp(cl/) 

since (again by (2.3)) the contribution of c > cq, cq a large positive constant, is 
negligible. If p is in the range indicated by 7^^r(a;), then 



^ ^ loga;p ^ ^ logx (i^q( Aog2^ 
logpr cL \ log a; 

log M = log2 X — log c — log L + O 

Therefore, for 2 < r < 7, 



logx 



r„(.)« E ^exp{-!^(l.0(j!gi:))x 



exp(cL— l)<p<exp(cL) 



X (|log2a;+ilog3a; + 0(l))} 
<a;exp|-cL(r-l)- fl + O 



2c \ \ log2 X 

so that the largest contribution is for r = 2. The function c+ l/(2c) attains its 
minimal value y/2 when c = 1/ hence 



%,r{x) < a;exp |-V2L ^1 + 



loggx 
logao; 



and we obtain 

(3.9) N{x) < a;exp |-v/21oga;log2a; (^1 + O } • 

This is somewhat weaker than the result implied by Theorem 1. But we saw 
that the main contribution to (3.7) comes from Tc^2{x), which is quite similar 
to the main sum appearing in the estimation (lower bound) of To{x) (cf. [12]). 
The only difference is that in the case of %^2{x) we have ^{xp~'^,2p), while in 
the case of Tq{x) the terms '^{xp~'^,p) appear. In the relevant range for p (i.e., 
ex.p{jL) < p < exp(2L)) this will not make much difference, so that using the 
arguments [12] we can obtain for N{x) the upper bound implied by the expression 
standing in (1.5)-(1.7). It is only for the sake of clearness of exposition that we 
gave all the details for the slightly weaker upper bound (3.9). This discussion 
completes the proof of Theorem 1. 
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4. Proof of Theorem 2 

For the proof of Theorem 2 we need an expression that is sharper than (3.5), 
which is only an upper bound. If n is counted by N(x) then 

a) Either P'^{n) divides n, and the contribution of such n is counted by Tq{x) 
(cf. (3.1)), or 

b) The number n is of the form n = pq^k,p= P{n), b>2, P{k) < p, P{k) i- q, 
where henceforth q (as well as p) will denote primes. Since n does not divide 
P(n)!, then (similarly as in Section 3) it follows that qb > P{n). Also, again as 
in Section 3, it is found that the contribution of p < exp(L/4) = Yi is negligible. 
The contribution of p > cxp(2L) = Y2 is also negligible. Namely in this case, since 
g^|n, q > P{n)/b 3> exp(2L)/ logo;, the contribution is clearly 

■Cx p"^ < a:exp(-|L). 

p>exp(2L)/ \ogx 

Therefore, if b) holds, we need only consider numbers n for which there is a number 
b >2 and a prime q e {p/b,p), such that n = pq'^k,p = P{n), P{k) < p, P{k) ^ q. 
The contribution of -P(/c) = q can be easily seen to be negligible, as can be also 
shown for the one pertaining to Pik) = p. Using the technique of Section 3 (or 
of [8]) it is seen that the main contribution comes from the sum with 6 = 2. If 
6=2, namely if n = pg^/c, p/2 < q < p and q^k with some prime p/2 < gi < p, 
then n is counted at least twice, but it is seen that the contribution of such n is 
^ a;exp(— |L). Therefore the contribution to N{x) of integers satisfying b) equals 

(4.1) 0(a;exp(-iL))+ J] J] * (£2'^^) • 

This is different from Ford [8, eq. (2)], who had W <7^ instead of the correct 

\E' (^^7 • This oversight does affect his result (1.1), which is not correct since the 
constant is not the right one (1 + log 2 should be replaced by 2), but the relation 
(1.3) remains valid, since both expressions with the ^f- function are of the same 
order of magnitude. It follows that our starting relation for the proof of Theorem 
2 takes the shape 

(4.2) iV(x) = To(x)+0(xexp(-|L))+ ^ ^ (^'^^) • 

By using the asymptotic formula (2.3) for \E'(x,2/), the prime number theorem in 
the standard form (see e.g., [14]) 

(4.3) 7r(a;) = y;i=/ -— + 0(a; exp(- ^1^)), 
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we write the sum in (4.2) as a Stieltjes integral and integrate by parts. We set for 
brevity 



R 



'\0g2X 



logx 

Then we see that the sum in question equals 



= (1 + 0{R))x r —p ( 



' 2.-1) d.dt, 



logt 



on making the substitution logy = .logt. Since both z and 1/z equal 1 + 0{R) 
in the relevant range, by the use of (2.5) we see that our sum becomes 

(l + 0(R))x P =^-1) / t-'dzdt 

Jy, tlogt \logt / Jl-log2/logt 



(4.4) 

= X(1 + 0(R)) 5 p 1^ 

Now we use (2.6) (with u — 3 replacing u, v = —2), setting u = logx/logt, 
Yi ^ t < to obtain that 

e2^("-i)p(« - 1) = e2^('^-i)p((« - 3) + 2) 
= e^^^^-^^-^^^'^-^Vlw - 3) ■ (1 + 0{R)). 

Since ^'(u) ~ 1/tt (m — > 00) we have 

g2^(«-l)-2^(«-3) ^ gO(l/n) = 1 + O = 1 + 0{R). 



Therefore the last integral in (4.4) equals 

(4.5) (i+o(fi))/%(^-3) 



dt 



log' t 



With the change of variable 

log X dt du 



u 



logt ' tXo^t log 



X 
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it follows that (4.5) becomes 

1 



{1 + 0{R)) y"''%(«-3)exp(^-i^^ du, 



S: = 

log a; 

1 / logx / logx 



log2 a; y log2 X 

But since up'{u) — —p{u — 1), integrating by parts we obtain 

S = ^(1 + 0{R)) r - 2) exp f-i^^ du, 

= (1 + 0{R)) P -{u - l)p'{u - 1) (^i^ exp d«, 
= (1 + 0(i.)) £ Pin - 1) (i^l^) exp (-i^) d., 
= loga;(l + 0{R)) j'J -p'{u) exp d«, 
= loga;(l + 0{R)) j'J p{u) {^^^ exp (-^) 

where we used several times that u x 1/R in the range of integration, so that 
lower order terms could be absorbed by the 0(i?)-term. Making again the change 
of variable log a;/ log t = tt, we obtain 



where (3.4) was used. This proves (1.8). 



X 
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5. Moments of the Smarandache function 

In this section we shall prove Theorem 3. We begin with the proof of (1.9). 
Recall that the Smarandache function S{n) denotes the smallest k E N such that 
n\k\. This implies that P{n) < S{n) < n, and if S{n) ^ P{n), then n does not 
divide P(n)!. Thus we may write 

^ S^(n) " ^ S^-ln) ^ ^ Win) 

2<n<x ^ ^ 2<n<x,S{n)=P{n) ^ ' 2<n<x,S{n)^P{n) ^ ' 

(5.1) ^ / \ 

2<n<a; ^ \2<n<x,n\^{n)\ ^ ' ) 

We have, by Ivic-Pomerance [12], 

(8-2) ^ 4? = "''^P {-^/2;L(x)(l + + o((!^)^ 

with L(a:) and Qrix) given by (1.6) and (1.7), respectively. On the other hand, by 
the argument that gives (3.5) and the proof of Theorem 1, 

^ P^(n) 

2<n<x,n\/P{ny. ^ ' 

« E E 

2<s<log x/ log 2 p<xi/s ^ 

< xexp |-(v'2(r + 1) + o(l))L(a;)| (x ^ oo). 

The proof of (1.9) follows then from (5.1)-(5.3), with the remark that a true 
asymptotic formula for the sum in (1.9) follows from [11, cq. (2.11)], since the 
major contribution to the sum in question comes from n for which S{n) = P{n). 
However, similarly to the comparison between (1.5) and (1.8), the advantage of 
(1.9) is that the right-hand side contains only elementary functions. 

The proof of (1.10) utilizes the following elementary lemmas, which generalize 
e.g.. Lemma 1 and Lemma 2 of De Koninck-Ivic [4]. 

Lemma 1. Let r > —1 be a fixed number, and let Cr,m = (m — l)!(r + I)"*" for 
m e N. For any fixed integer M > 1 we have 

^/ \logx log^a; Vlog^+'x 
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Proof. Follows by partial summation from the prime number theorem in the 
form (4.3). 

Lemma 2. For a given fixed real number r > and given natural numbers j 
and m, we have 



(5.5) 



1 



n< 



J]C(^)(r + l) 



fc=0 



-3 



1 



+ 



k J {\ogxy+^ yiog^+'^+^x 



1 



Proof. Set I = l{x) = exp(vloga-). Then 



y I = y L_ 



log n 
logx 



+ 0(0- 



The assertion of the lemma follows if we use the binomial expansion 
\ogn\ _v-^_pfe 
fc=o 



1 



log X J V ^ / V ^ / \ V ^ 

and (C(s) denotes the Riemann zeta-function) 

oo 

log^ n • n"'-^ = (-1)'=C^''^ + 1) (r > 0; /c = 0, 1, . . . ). 



Now we turn to the proof of (1.10). We have, by Theorem 1, 
(5.6) 

2<n<a; 2<n<a;,S(n)=P(n) 2<n<K,S(n)5!^P(n) 

Y p»+o(x'- 5; 1 

2<n<a;,S(n)=P(n) \ 2<n<x,n\P{n)\ 

= Y P''{n) + o\x'-+''e-K^{-{y/2 + o{l))L{x))^ {x ^ 00). 

2<n<x 

Observing that ^{x/p,p) = [x/p] if ^/x < p < x, it follows that (A;, n e N, r > 0) 

Y p'^in)^ Y p'-T.p'^(l^p) 



2<n<a; 



pk<x,P{k)<p p<x 



(5.7) 



p<\fx 



^<p<x 



o(xY p"~') + Y p' 

p<\fx ^/x,<p<x 

o(xi+-/2) + Y p"- 



pn<x 
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By using Lemma 1 and Lemma 2 we infer that 



with 



(5.8) 



Ef"" E E 

pn<x n<ix p<x/n 



n<^x 



+ 



X 



J J 



r+l 



X 



-E-,Ec'-(..i)(7)p^.o 



X 



r+l 



log-^^^ X 



fc-i 

ak,r := J2^k - i - l)!(r + l)^-'=C^^Hr + 1) 

i=0 



— I 

k — i 



The assertion (1.10) of Theorem 3 follows then from (5.6)-(5.8), clearly one has 
ai,r = C(r + l)/(r + l). 



(5.9) 



Finally we remark that, for r e R fixed and a; — > 00, 

S{n) 



E (^X =x + o{xeM-iV2 + omm}- 



2<n<x 



Namely, if n = p^^ . . .p^^ with pi = P{n) is the canonical decomposition of 



then n divides 



{(ai + a2 + . . . + ak)pi}\ < \ P{n) 



+ 1 ! 



logn 
1^ 



Thus we have P{n) < S{n) <^ P{n) logn, and one obtains without diflficulty (5.9) 
by the method of proof of Theorem 3. 
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